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ABSTRACT 

Modified Newtonian Dynamics (MOND) can fit a broad range of galaxy kinematic 
data, but struggles with clusters of galaxies. MONDian clusters need dark matter, and 
here we test the 11 eV/c^ sterile neutrino - used to fit the first three acoustic peaks 
of the cosmic microwave background - by investigating their equilibrium distributions 
in 30 groups and clusters over a wide range of temperatures. We do this by first 
taking the known sterile neutrino density, necessary for hydrostatic equilibrium of the 
intracluster medium (or to produce the observed lensing map). Then, we solve for 
the sterile neutrino velocity dispersion, needed for their own hydrostatic equilibrium, 
through the equation of state for a partially degenerate neutrino gas. The velocity 
dispersion is a unique, continuous function of radius determined by the density and 
mass of the sterile neutrinos particles. Knowing both the sterile neutrino density and 
velocity dispersion tells us the Tremaine-Gunn phase-space limit at all radii. We find 
that all 30 systems serendipitously reach the Tremaine-Gunn limit by the centre, 
which means a portion of the dynamical mass must always be covered by the brightest 
cluster galaxy. Interestingly, the typical fitted K-band mass-to-light ratio is unity and 
at most 1.2, which is very consistent - although leaving no margin for error - with 
stellar population synthesis models. Amidst the sample there are several special cases 
including the Coma cluster (for which dark matter was first proposed), NGC 720 
(where geometrical evidence for dark matter was found) and the bullet cluster (where 
dark matter - of some kind - in clusters was directly proven to exist). We demonstrate 
that 11 eV jc? sterile neutrinos are unlikely to influence spiral galaxy rotation curves, as 
they don't influence even some very massive early-types (NGC 4125 and NGC 6482). 
Finally, we conclude that it is intriguing that the minimum mass of sterile neutrino 
particle that can match the cosmic microwave background is identical to the minimum 
mass found here to be consistent with equilibrium conflgurations of MONDian clusters 
of galaxies. 



1 INTRODUCTION 



Sterile neutrinos are hypothetical additions to the standard 
model of particle physics. They are right handed, neutral 
leptons which interact only via gravity, which earns them 
the "sterile" prefix, contrary to the active neutrinos which 
also participate in the weak interaction. Gravity aside, ster- 
ile neutrinos also can interact with the active neutrinos via 
the quantum mechanical phenomenon of neutrino oscilla- 
tions. This behaviour has bee n investigated by th e Liquid 
Scintillator Neutrino Detector ('Aguilar et al."200l') and the 
MiniBoone experiment (,MaI toni & Schwetz 2007i ) . However, 
no concrete evidence was convincingly found to suggest the 
existence, or non-existence, of sterile neutrinos from the 
disappearance of active neutrinos. The basis of some ap- 
peals to sterile neutrinos are for aesthetic reasons, since the 



active neutrinos are entirely left-handed; and others make 
use of them in the so called "see-saw mechanism" which 
can give rise t o the small masses of the active neutrinos 
([Lindner et al.l i2002). Here we make no claims of a deeper 
theory for sterile neutrinos, but rather continue to investi- 
gate a startling coincidence. 

A rece nt anal y sis of the Cosmic Microwave Background 
(CMB) by I Angus! (|2009l ') has demonstrated that the acous- 
tic peak s in the angular po wer spectrum as measured by 
WMA P l|Dunkley et al.l2009l ) and ACBAR l|Reichardt et all 
120091 ) can convincingly be generated by a single, thermal (by 
virtue of neutrino oscillations in the early Universe) sterile 



neutrino (SN) with mass 



= 11 eVI(?. This SN is a 



straight substitution for the cold dark matter (CDM) of the 
concordance cosmological model (jSpergel et al.ll2007l ). such 
that = 0.0205m^^ = 0.225, f^b = 0.047 and the spec- 
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tral index Us — 0.965. However, whereas the unknown CDM 
particle can condense and form structures on virtually any 
scale, which puts it at odds with certain observations - like 
the lack of substruc ture as compared t o the cosmoloRical 
N-body simulation s (|Klvpin et al.lll999l : iMoore et al.lll999l : 
iKroupa et al.ll2005ll and the e xistence of DM core s in galaxy 
halos l|Gnedin fc Zhaoll2002l : iGentile et all 120041 ') - the hot 
SN free streams out of galaxies. This makes it a non-starter 
as a traditional dark matter candidate, since it would render 
spiral galaxies (like the Milky Way) and dwarf galaxies (like 
Draco and Ursa Minor) incapable of retaining their stars. 

Fortunately, there is a well studied alternative to the 
standard lore of gravi ty called Modified Newtonian Dy- 
nami c s (or MOND; see iMilgromll 19831: ISanders fc McGaughl 



l2002l : iBekensteinI l2006l : iMilgromI l2008bl " for reviews) which 
has an uncanny knack of predicting the dynamics of galaxies 
from a simple relation using only the baryonic matter den- 
sity. The paradigm of MOND introduces a new physical con- 
stant with dimensions of acceleration, Qq = 1.2 x 10~^''ms~^, 
around and below which dynamics do not follow from stan- 
dard Newtonian theory. In particular, the true modulus of 
gravity, g{r) = T4(r)^r~^, is not linearly related to the New- 
tonian gravity, gn{r) = GM{r)r~^ , but instead 



9 = ^9n 



1 + Wl 



4ao 



(1) 



This ensures adherence to two critical axioms: that gravity 
is Newtonian in regions of strong gravity, and that when 
g « Oo, g oc 1/r meaning rotation curves are fiat at the 
periphery of spiral galaxies. 

The additional gravity afforded by M OND, replaces 
the need for D M in dwarf spheroidals (lAngud 
Milgroml 1 19951: ISanchez -Salcedo fc Hernanded 



Serraet al. '2009*), spiral (e.s. 

19981 : ISanders fc Noordermeej 12007 
200, 



2008 
2007 



McGaugh fc dc Blokl 

. _ . . : [Famaey & Binney 

McGaughl l2008l) and X-ray dim ellipti cal galaxies 



ijMilgrom fc Sanders! l2003l : lAngus et al ] l2008bl ) often with 



remarkable accuracy. Furthermore, rotatio n curves of tidal 
dwarf galaxies have been observed by iBou rnaud et al 
2007) an d the indepen dent analysis by [Gentile et al 



2007^ and lMilgronJ (|2007l ') reveal not only their consistency 



with MOND (with zero adjusted parameters), but also 
constitutes a direct falsification of dark matter being made 
only of CDM i n galaxies: assuming the data are reliable. 
What is more, iKroupa et al.l (|2005r ) have advocated that 
the Milky Way's dwarf spheroidal galaxies are also tidal 
dwarf galaxies. This convincing demonstration of MOND's 
predictive ability does not, however, extend to the realms 
of clusters of galaxies. 

An obvious difference between galaxies and clusters of 
galaxies is simply the scale involved. If a typical galaxy is 
one or two tens of kpc across, clusters have accurate mea- 
surements of the gravitational potential (through the intr- 
acluster medium or weak gravitational lensing) on scales 
ten times that and therefore volumes many thousand times 
greater. Studies of the dynainics o f groups and c l usters 
in MOND (e . g. |T he fc White 'l988'; Sanders 'l99 4|, Il999l: 
Aguirre et al] 12001.: Sanders ,2003; Poin tccoutcau fc Silkl 



20051: lAngus et al.ri2007l : [SandersI I2OO7I : iFefac et al.1 [2OO8I : 
Milgroml l2o68bF have all shown that there is a huge cen- 
tral mass deficit and lAngus et al] (|2008al ) demonstrated 
that active neutrinos, even at the experimental maximum 



(~ 2 eV/c^), cannot clump densely enough. A more plausi - 
ble solution is the 11 eV/c^ SN proposed bv I Angu j (|2009l 'l 
to fit the first three peaks of the CMB. 

The feasible particle mass of SN is fully determined by 
fitting the CMB power spectrum, and there is virtually no 
freedom available above 5% of the 11 eV/c^ mass. It is the 
particle mass, and particle mass alone, that sets the prop- 
erties of the SNs in clusters and galaxies. Therefore, it is a 
highly intriguing corollary that the mass required to match 
the CMB is in the tiny range of neutrino masses (perhaps 
11 — 12 eV/c^) that can both free stream out of galaxies 
(MOND does not require any DM in galaxies) and clump 
densely enough in galaxy clusters to account for the serious 
mass deficit exposed there. 

We now have the basis of a predictive cosmological 
model, where we have made two positive trades. Firstly we 
exchange Newtonian dynamics for Modified Newtonian Dy- 
namics, which helps explain in detail the origin of the mass 
discrepancy in all galaxies, and with fewer freedo ms than the 
often contrived CDM halos (as d emonstrated bv iMcGaughl 
l2005l : lKuzio de Narav et aLlbOOgl ). Secondly, we swap CDM 
for an 11 eV/c^ SN. The added bonus of SNs over CDM, is 
that knowing only the particle mass gives fixed predictions 
for the CMB and for structure formation (with a small de- 
pendence on the /i- function in only the latter case). 

In the very early Universe, neutrino decoupling (both 
active and sterile) occurs at a temperature of kT ~ 1 MeV. 
This means that the 11 eV/c^ SNs (as well as the sub 
eV/c^ active neutrinos) are ultra-relativistic during decou- 
pling, which freezes in their Fermi-Dirac distribution an d 
their cosmological abundance is fixed (e.g. |Peacoclj|l999l ). 
Typical CDM candi dates, like 100 GeV neutralinos (e.g. 
iHofmaim et al.ll200ll) or > 100 keV/c^ SNs (see the review 
bv lBovarskv et al.ll2009l ) are non-relativistic during decou- 
pling, and self-annihilations must be used to tune the cos- 
mological abundance. This applies to all CDM candidates 
which means knowing the mass of a CDM particle tells us 
the mass and not the cosmological abundance. Free param- 
eters, like the interaction cross-section, must explain why 
they have the correct cosmological abundance. 

In addition, while searching for an explanation of the 
anomalous low energy excess of e lectron neutrinos o b served 
by the Miniboone experiment, iGiunti fc Lavedeil (|2008l ) 
found agreement with the data by postulating a perfectly 
plausible renormalisation of the original flux of muon- 
neutrinos and oscillations (that are energy dependent) from 
electron neutrinos to 11 ± 7 eV/c^ SNs. This hypothesis 
clearly warrants further investigatio n, and will be t estable 
with the T2K neutrino experiment (|Hastingsl |2009| ) which 
will have a near detector at L = 280 m giving it a similar 
L/E (E being neutrino energy) as MiniBoone (M. Laveder, 
private communication) . 

In this paper, we seek to investigate the 11 eV/c^ SN's 
influence in clusters of galaxies in MOND. Primarily, we 
need to ascertain that every cluster has an equilibrium dis- 
tribution of neutrinos that can account for all the gravi- 
tating mass. This is by no means guaranteed since there 
is a maximum density set by the Tremaine-Gunn limit 
llTremaine fc Gunnlll979l : see ^3.4^ depending on the veloc- 
ity dispersion of neutrinos and the mass of an individual 
particle. If the required density of the equilibrium models 
exceed the Tremaine-Gunn limit (or if the required AI/L 
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of the BCG to keep the SN density below the Tremaine- 
Gunn hmit is too high), they would be excluded. To this 
end we have taken a relatively large sample of 30 relaxed 
galaxy groups and clusters to gauge their consistency with 
the missing mass problem of clusters of galaxies in MOND. 



2 DATA 

We take the temperature and density profiles of the intra- 
cluster medium (ICM) in 7 clusters of galaxies with tem- 
peratures in th e rang e of 1.8 — 9 keV from the sample of 
IVikhhnin et~a r ('2006'). a further 13 _groups of galaxies from 
the s ample of G astal dello et all (|2007l ). the cluster A 2589 
from IZappacosta et alT i 20061 ) and 3 X-ra y bright, isolated 



early type galaxies from the sample of iHumphrev et al 



This covers the full sample used in 



Angus et al 



2008al ). but in addition, we make an analysis of the two 



clust e rs comprising the bullet cluster (iClowe et"al ] l2004 
l2006l : iBradac et alTl2006l : lAngus et al] |2007| ) working from 
the NFW halos fitted to the weak-lensing convergence map. 
We also take a fitted NFW profile f or the cluster A 168 9 
from the strong lensing analysis of iHalkola et al.l (|2006l ) , 
take 3 estimates for the density profile of the Coma clus- 
ter (jGavazzi et al ] |2009l: iKubo et al.l | 2007l 'l. a pair of clus- 
ters in the Lynx field fjee et al.''2006) and the fossil group 
RXJ 1416 jKhos roshahi ct al. 2006 ). 

For the lGastaldello et al.l (|2007l ) sample, we only had X- 
ray data for 10 of the 13, so for the other 3 (A 2717, IC 1860, 
MS 1 160) we started from their fitted NFW (|Navarro et all 
1 19971 ) profiles and included t he brightest cluster galaxy. 

In lVikhlinin et aTl ()2006l ). the authors give fully analyt- 
ical descriptions of the ICM temperature and density allow- 
ing us to solve the equation of hydrostatic equilibrium to 
give the dynamical mass of each cluster (see Similarly 
for the groups of galaxies, we were provided with (D. Buote 
private communication) high resolution data that gives the 
ICM density and temperature as functions of radius. We fit- 
ted /3-models to the ICM de nsity and have fi tted the same 
analytical models defined in IVikhlinin et al.l ({2006 ) to the 
temperature. These models are rather sophistocated, many 
parameter models that maximise the accuracy of the fits. 

The masses of the brightest cluster galaxies 

i BCG s) are taken ei t her f ro m the Gasta ldello et al.l 
20071): IVikhlinin et all (|2006l ): iHumphrev et al. ( 2006); 
IZappacosta et ahl ( 20061) K-ban d lum inosity estimates 



the observations of iLin fc Mohil (|2004l ) or f rom K-band 
magnit udes on the NED database as per lAngus et al.l 
(|2008al ). 

Naturally, there are errors associated with the best fit 
density and temperature profiles, which we do not explicitly 
investigate here. It is enough to say that the observational 
uncertainties surrounding the Newtonian dynamical mass 
(from which everything is deduced) are no larger than the 
uncertainties in the triaxiality of the cluster or the interpo- 
lating function used to find the MONDian dynamical mass, 
so we waste little time debating them. Moreover, we use a 
sample containing 30 of some of the most relaxed systems 
so that a general consensus can be reached. 



3 SOLVING FOR STERILE NEUTRINO 
EQUILIBRIUM CONFIGURATIONS 

3.1 Intracluster medium hydrostatic equilibrium 

The density of 11 eV/c^ SNs (or any DM candidate) in 
MONDian clusters of galaxies required to provide hydro- 
static equilibrium of the ICM is given by the following steps, 
as per lAngus et al] (|2008al ). Firstly, we take the observed 
density and temperature of the ICM, px{r) and kTx{r) re- 
spectively, and numerically take their logarithmic derivative 
with respect to radius to find the true gravity as a function 
of radius 



air) 



-kT^r) 



dlnp^r) _^ d\nkT4r) 



(2) 



d In r d In r 

where w = 0.62 is the mean molecular weight of the ICM. 
The gravity simply is related to the total MOND enclosed 
mass by 

M^{r)=^G-'g{r)p{g/a,), (3) 
and the interpolating function is the simple one 



Ka/ao) 



(4) 



1+ g/tto' 

(see iFamaev et al ] |2007l for a discussion of how it fares in 
spirals). This single line is the only stage at which MOND 
is involved. 

After subtracting the mass of the ICM {Mx{r) = 
Attv^ Px{f)df), we are left with the SN mass distribu- 
tion and the unsubtracted BCG, Mtcg + M^^ (r) = Mm(r) — 
l.l^Mx(r), where the mass of the brightest cluster galaxy 
{Mbcg = M/Lk X LK,bcg)- The 1.15, which is of virtually no 
consequence, multiplies the ICM mass to include the contri- 
bution of galaxies in the cluster. We invert this to give the 
SN density, there we ignore the BCG until later, although 
it is vitally important, for reasons that will become obvious 
later. 



dr 



3.2 Sterile neutrino hydrostatic equilibrium 



(5) 



At this point we have a deduced density of SNs in a cluster 
of given temperature profile: at least it is the density exactly 
required for hydrostatic equilibrium of the ICM (or in the 
case of lensing, to create the observed convergence map). In 
addition to this constraint, there is an equation for hydro- 
static equilibrium of the SNs themselves, which is crucial to 
this analysis. This equation invokes the equation of state of 
neutrinos (active or sterile; see e.g. Sanders 2007) and what 
this gives us are equations that define the density and pres- 
sure of a partially degenerate neutrino gas, coupled via the 
hydrostatic equilibrium relation 

^P..=-p.Ar)g{r). (6) 

To express the density, we start with the equilibrium 
occupation number 

1 



/ = 



1' 



(7) 



exp {x - x) 

where x and x are respectively the ratios of neutrino energy 
and chemical potential to temperature. A large positive or 
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negative x denotes strong degeneracy or non-degeneracy re- 
spectively. From this we find the phase space density 



(8) 



Here, h = 4.136 x 10~^^ey.s is Planck's constant and gvg=2 
takes into account the anti-particles, whilst remembering 
that neutrinos have only a solitary helicity state. One im- 
mediately sees that there is an absolute upper limit to the 
allowed density of SNs in phase-space, corresponding to full 
degeneracy (x = +oo) and known as the Pauli limit. The 
starting momentum distribution of neutrinos in the early 
universe corresponds to half this limit. 

It should be clear here that in groups and clusters of 
galaxies, SNs of 11 eV/c^ mass are non-relativistic since they 
are travelling at velocities of the order 100-1000 km s^^ (e.g. 
Fig H central panels). Therefore, neutrino energy is related 
to momentum by _E = /2mv^ . The number density of neu- 
trinos is given by J fd^p and multiplying this by neutrino 
mass (m^j, =11 eV/c^), converting momentum to energy 
and temperature to velocity dispersion (fcT^^ — rriv^fJ^^) we 
find our equation for neutrino mass density 



(9) 



where the SN velocity dispersion and temperature are ct^^ 
and kTv, respectively and -Fi/2(x) = x^^^ fixjdx. In a 
non-relativistic neutrino gas, the pressure is equal to 2/3 the 
internal energy per unit volume (Uv^ = J E{p)fd"*p) giving 



(10) 



where F3/2(x) = x^^^ f{x)dx. 

It emerges that there are two variables here that must 
be set in order for the neutrinos to exist in hydrostatic equi- 
librium. Primarily, x('")i niust be fixed so that the fixed den- 
sity of Eq[5]is matched by Eq[9l In case this is not obvious, 
we solve Eq|5]for -Fi/2(x)) for which there is a unique, con- 
tinuous x('") once the neutrino velocity dispersion is given. 
There is no independent method, apart from cosmological 
simulations of the collapse of the baryonic plus SN two-fluid 
mixture to estimate the chemical potential, X) so we choose 
here to fit it to the data. Whether detailed numerical sim- 
ulations of cluster formation in a MOND cosmology will 
reproduce this chemical potential will be a crucial test of 
the equilibrium models presented hereafter. 

This x('') is then transfered to Eq 1101 and although 
x(r) is used to balance p^(r) and ensure it remains un- 
changed when is varied, it infiuences the pressure, P^^, 
non-linearly. Secondly, Eq [S] must be satisfied, however, it 
is not satisfied if we make the assumption that the sterile 
neutrino velocity dispersion (VD) is identical to the ICM 
VD. In general, there would be too much pressure because 
the neutrinos would be too hot for their given distribution. 

A priori, there is no reason that the 11 eV/c^ SNs 
should have precisely the same temperature as the ICM. 
Both fluids are orbiting in the same gravitational potential, 
but there is no rapid exchange mechanism to bring them into 
mutual equilibrium, as there was in the very early Universe. 
In fact, this is likely to be related to the formation epoch 
of the cluster, where the SN halos would have formed at 
relatively high redshift while the Universe was more dense. 
Later, the ICM would have fallen into the deep potential 
well, created by the SN halo, and therefore it is logical that 



the ICM should be hotter than the SNs. In addition to prob- 
ing the x('")'s fitted, numerical simulations of galaxy clusters 
in MOND will, in time, be able to tell us if the difference in 
neutrino and ICM VD is realistic. 

From Eqs[9]&[T0]one can see that the SN VD influences 
the product of the density and gravity (the right hand side 
of Eq |6} non-linearly with respect to the gradient of the 
neutrino pressure. Essentially, to force the SN distribution 
into hydrostatic equilibrium, with its fixed density, we can 
modify the SN VD. However, as we shall discover in the next 
section, this cannot be achieved by a simple scaling from the 
ICM VD, but instead by solving for the neutrino VD, cr^^, 
as a unique, continuous function at every radius. 

3.3 Solving for a unique sterile neutrino velocity 
dispersion 

In principle, for hydrostatic equilibrium of the SNs, we need 

US 



ar 



1 and writing this out after some algebra gives 
<(0 3 Fi/2(r) 



5F3/2(r-)f F3/2(r) W F.^^ir) 



5W.(11) 



Using the substitution e(r) = al^{r), leading to e'(r) = 
2a,y^{r)al^{r) and the integrating factor -F3/2('')'''*, we can 



reduce this first order linear differential equation to 



(^I'sii^) = g-F3/2(r)" 



Fi/2{f) 



3/2 



(12) 



To find pressure and density in the first instance, we 
must submit a trial to Eq 1121 typically we try ^o-^, 
where the ICM VD is defined by 

1/2 



(13) 



Here, c = 3 x 10^ kms"^ is the speed of light, kT^ is the ICM 
temperature in keV, w — 0.62 is again the mean molecular 
weight and rrip — 9.38 x 10^ keV is the mass energy of a 
proton. From the trial solution, we iterate rapidly towards 
convergence. The final a^^ (r) gives hydrostatic equilibrium 
to better than 1 per cent at all radii and most importantly is 
unique - set only by the SN mass and its indirectly observed 
density from Eq[Sl 

3.4 Phase space constraints 

We now have a unique correlation between the density of 
SNs and their velocity dispersions. This is important because 
these two variables define the phase space distribution of the 
SNs, to which there is a fundamental limit. 

Liouville's theorem states that (in the absence of en- 
counters) flow in phase space is incompressible and that each 
element of phase-space density is conserved along the flow 
lines. However, this only applies to the fine-grained phase 
space density of an infinitesimal region. Rather, for the ob- 
servable, which is the coarse-grained (macroscopic) phase- 
space density, we simply must not exceed the maximum of 
the fine-grained one. 

Thus, the SN phase-space density must not increase 
during collapse, from its starting value of ^g^h~^mf,_, (which 
is half the Pauli degeneracy limit), to its current maximum 
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value of pi/j (r)[27rCT^^ (r)]^^'^ (where we assume the veloc- 
ity distribution is locally Gaussian everywhere with dis- 
persion a^,). This limit is called the Tremaine-Gunn limit 
IITremaine fc Gunnlll979l : hereafter TG) and rearranging it 
in terms of the critical density for an 11 eV jc? SN, where 
1 eVjc^ = 8.9x 1O"^^M0 and s^c^ = 9.18x lO'^'^pc^ (where 
s is obviously 1 second) gives us 



p..,TG(r)= 3.15 xlO«(^^y 



Mqvc 



(14) 



The important thing to bear in mind here is that this is 
the exact value of the TG limit because it is calculated from 
the unique, derived SN VD exactly necessary for hydrostatic 
equilibrium and not by assuming some relation between the 
ICM and the neutrinos, as is often the case ( Angus et al.l 
l2008al : iNataraian fc Zhadl20o4 iGentile et al]|2008l ). There- 
fore, we have the TG limit as a function of radius, since we 
know the exact value of the VD from solving Eg 1121 

Let us note however that departures from Gaussianity 
in the velocity distribution could perhaps somewhat affect 
this physical upper limit on the SN density. 

The consequence of the TG limit is that, therefore if the 
density of SNs required for cohesion of the cluster gas ex- 
ceeds TG limit (in particular at the centre) then the MOND 
plus 11 eV ji? SN hypothesis would be ruled out. 



4 RESULTS 



In Figs 11161 we present the densities, VDs and enclosed mass 
profiles for the 30 groups and clusters. In the left hand pan- 
els we plot the SN density (solid line type), ICM density 
(dotted) and the TG limit of the cluster (dashed) against 
radius. In the middle panels we plot both the observed ICM 
(dotted) and the derived (from Eq[T2]) SN (solid) VD as func- 
tions of radius. In the right hand panels, we plot each clus- 
ter's enclosed Newtonian dynamical mass (dashed), MOND 
dynamical mass (dot-dashed), BCG total mass with unity 
MjL (sohd) and ICM mass (dotted). 

There are a few salient features to observe: primarily, 
the SN VD (or temperature) is in most cases 20-50 per cent 
lower than the ICM VD. The simple explanation, alluded to 
earlier, is that the SN halos presumably formed at relatively 
high redshift, through only their mutual gravitation. On the 
other hand, the ICM fell from large distances, through the 
already present, deep potential well of the SNs and thus 
had greater potential energy to transfer to kinetic energy, 
although this still has to be demonstrated with numerical 
simulations. 

Taking a closer look at the SN densities and VDs, they 
have conspicuous kinks around 20-30 kpc. This radius, rtg, 
is where the density reaches the TG limit, . If no phase space 
limit existed, the SN density would continue to increase to- 
wards the Pauli limit and the neutrino equation of state 
would begin to substantially change from non-degeneracy 
(large negative x) Pvs oc Pv^'^ts degenerate one 

(large positive x) Pvs c< pV^ (where the TG limit occurs 
at X ~ 0.35). However, since we must adhere to the TG 
limit, we impose it by making the density equal to the TG 
limit. The density and VD both increase towards the centre 
as they must conform to the duel constraints of remaining 



at the TG limit and satisfying hydrostatic equilibrium. Al- 
ternatively, if we had ignored the TG limit, the VD would 
have crashed towards the centre of the cluster; as less classi- 
cal pressure would be required to prop up the SN halo. This 
is simply because gradually more Fermi pressure would be 
available. It is this behaviour that causes the kink since we 
force the density to be equal to the TG limit, the required 
VD for hydrostatic equilibrium must rise, whereas it was 
falling prior to this. 

Since the SN halo of every group and cluster reaches 
the TG limit at the centre when there is still a discrepancy 
between the central SN density and the required total den- 
sity for hydrostatic equilibrium of the gas, the BCG plays 
a crucial role in determining how satisfactory our 11 eV/(? 
SN hypothesis is. If no BCG existed in any of our clusters, 
they would immediately fail. 

We know the luminosity of most BCGs in the K- 
band which gives an excellent indication of mass (see e.g. 
i Bell fc de Jondl200ll: [Conrov et al.ll2009l : iGastaldello et al.l 
'2007': Hu mphrev et al.l |2006| ) and we know a few in other 
bands (B,V,R) or not at all. Therefore, we must modify the 
M/ Lk (given in table 1) to exactly match Mm at rtg, which 
basically means the BCG is picking up all the "slack" left 
when the SN reaches the TG limit and can no longer account 
for the full dynamical mass. 

There are some points to bear in mind, firstly the high- 
est M /Lk used was 1.2 and the lowest 0.1. In the latter case, 
of a very low M / Lk , this can be due to the total luminosity 
not being enclosed by rtg which is often ~ 20 — 30 kpc). 
In the case of M/ Lk ~ 1, this can be considered a fitting 
parameter, since there is no freedom in the cluster to have 
a M/Lk lower than the value used because the SN density 
reaches the TG limit at the centre of every cluster. Further- 
more, the M/Lk cannot be significantly larger or it would 
be in disagreemen t with the typical M/Lk demonstrated 
bv lBeU fc de Jongj l|200ll '). In the cases where M/Lk < 0.8 
it need not be considered a fitting parameter, since there is 
freedom for the M/Lk to be larger and simply for the SN 
to be lower. Nevertheless, the TG limit will still be reached 
in every cluster, the only difference will be that if we are un- 
derestimating the BCG mass then rtg will simply be lower 
for that particular cluster. The only reason we fix the den- 
sity to be equal to the TG limit for all radii smaller than 
rtg is to highlight the maximum amount of luminosity the 
BCG could lose (e.g. if observations are incorrect) and still 
provide the central density required. 

Generally, it is interesting to note that the SN halos 
needed to fit galajcy clusters here have a density slope sim- 
ilar to that of the ICM in the central parts, but becoming 
sharper at intermediate distance, which is accompanied by 
a relatively flat velocity dispersion. At the edges of the clus- 
ters, the ICM density becomes larger than the SN density 
(which falls to zero) and there is an apparent sharp decrease 
of the sterile neutrino velocity dispersion to zero also. This 
is merely a numerical artefact of the sterile neutrino density 
being set to zero at the edges of clusters, under the assump- 
tion of hydrostatic equilibrium. Since the sterile neutrino 
density will in reality fall to a very small number, but still 
greater than zero, the velocity dispersion will actually be 
isothermal, as we would expect. 
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4.1 Individual systems 

Below we discuss some pertinent observations about indi- 
vidual groups or clusters. 

4.1.1 The bullet cluster 

Analysis of the CMB strongly favours the hypothesis that 
non- baryonic DM exist s in the Universe, and the b ullet clus- 
ter |Clowe et all l2006h compounds it (although iMilgronil 
l2008tJ has suggested that the DM of clusters could be ultra 
cold and coUisionless clumps of molecular gas which would 
satisfy the constraint of this particular cluster) . As two giant 
galaxy clusters crashed into each other at incredible speed 
iMarkevitch et alj |2002| . but see also lAngus fc McGaughl 
I2OO8I for why this might pose a problem for AC DM), mu- 
tual ram pressure from the 2 ICMs dragged one another 
out of the clusters leaving only the galaxies and DM to 
pass through and emerge on the opposite sides on the sky. 
A weak-lensing reconstruction required two DM halos to 
overlay the positions of the galaxies with NFW parame- 
ters for the main and sub cluster respectively - M200 = 
15.0,1.5 X 1O"M0, ''200 = 2.1, 1.0 Mpc and concentrations 
c= 1.94, 7.12. 

In this case, where we are directly given the Newtonian 
mass profile, we begin the procedure at EqOand follow the 
same steps. As with the other clusters, in Fig |6] we show 
that these two clusters both reach the TG limit smoothly 
at the centre. However, notice we must add a 3.5 x 10^^ A/0 
stellar mass to the sub cluster (labeled "Bullet 2") to cover 
the dynamical mass in the central 30 kpc, which i s very 
significant. The stellar mass quoted in IClowe et all (|2006l ) 
at the weak-lensing peak of the sub cluster is (5.8 ± 0.9) x 
10^^ Mq within 100 kpc, from I-band observations assuming 
a M /Li of 2. We need only add a trivial stellar mass of 1.0 x 
10" M0 to the main cluster, which is interesti ng since one 
can se e from the left hand panel of Fig 1 from lClowe et all 
l|2006t ) that there is no obvious BCG candidate associated 
with it. In fact, the two giant ellip ticals of the main cluster 
highhghted by IClowe et all (|2006l ) are roughly 50 kpc (to 
the northern one) and 75 kpc (to the eastern one) from the 
lensing peak, but can easily offer the necessary stellar mass. 
On the other hand, the centre of the BCG of the sub cluster 
is only roughly 25 kpc from the lensing peak and significant 
light is spilling over within 10 kpc. Furthermore, since BCGs 
usually dominate the stellar mass in the central 100 kpc of 
clusters, the majority of the 5.8 ± 0.9 x 1O^'^M0 is likely 
associated with it, making the 3.5x 1O^'^M0 plausible. Notice 
also that the 3.5 x 10^^ is not required within 10 kpc, 
but rather by 30 kpc. Within 10 kpc less than 1.5 x 10^^ M0 
is required. 

4.1.2 RGH 80 

Another intriguing point about BCG masses is the one 
used for the group RGH 80, which has two equally mas- 
sive central galaxies: NGC 5098a and NGC 5098b with 
Lk = 2.9 X 10" 1/0 and Lk = 2.4 x IO^ Lq respectively. In 
our previous paper l|Angus et al.ll2008al ) we were only inter- 
ested in the central 100 kpc, therefore, were able to combine 
the two galaxy luminosities together. Howe ver, only one of 
the galaxies is at the very centre (see Fig 9 of lMahdavi et all 



I2OO5I or Fig 1 of iRandaU et al.ll2009l ). and the other has a 
projected separation of around 50 kpc and possibly a con- 
siderable line of sight distance. Before removing the second 
galaxy, the TG limit was considerably larger than the max- 
imum central density, but after discounting it, the TG limit 
is reached at 20 kpc. 

4.1.3 Clusters with no BCC luminosity: A 1689 and 
A 2390 

We found no galaxy luminosities for the two clusters A 1689 
and A 2390, but they require BCGs with total masses of at 
least 2.0 and 18.5 x 10^^ Mq respectively, which is a predic- 
tion of this model. Interestingly, A 1689 has been used ex- 
tensively to argue against CDM by Broadhurst & Barkai3 
(|200^ because the observed NFW concentration parameter 
is considerably larger than that expected from cosmological 
simulations. We find the equilibrium model with 11 eV/c? 
SNs nicely reaches the TG limit at 20 kpc, even though it 
has one of the highest central SN densities of all our sample. 

4.1.4 Lynx 

The two clusters in the Lynx field are in the process of 
merging (like the bullet cluster) and thus the mass profiles 
are potentially overlapping. Nevertheless, their offset seems 
to be sufficiently large to get a reasonable estimate from 
weak lensing and th is has also been sanity checked with X- 
ray hydrodynamics (|jee et al.ll2006l) . Both the larger cluster 
and smaller cluster are fitted with the same NFW profile: 
M200 = 2.0 X 10^* M0, ^200 ~ 0.75 Mpc and concentration 
c = 4. The galaxy luminosities within 500 kpc of each of the 
two lensing peaks are 1.5 and O.8xlO^^L0 in the B-band. 
We require 6.0 x lO^^Af0 for each BCG, which is easy af- 
fordable by the luminosity of both clusters because typical 
B-band M/ Ls can range between 5 and 10. 

4.1.5 Coma 

The historical significance of the Coma cluster with respect 
to the dark matter problem is probably far more signifi- 
cant than its scientific significance in this present case, but 
we include it out of curi osity. It was originally analysed in 
MOND bv lThe fc Whitel (19M) and was actually concluded 
to be more or less consistent with MOND, although the 
galaxies required radially biased orbits and the accelera- 
tion constant, ao had to be increased by at least a factor 
of 2 from the one used here. The problem with Coma is 
that it is not particularly r elaxed and also it s sph ericity is 
questionable. For example, iNeumaim et all (|2003l ) showed 
that there is ongoing merging, which makes measurements 
of the central mass profile uncertain. Therefore, the only 
way to get a decent estimate of the dark halo of Coma 
is to use weak lensing (like we have done with the bullet 
cluster) , but even then the Eissumption of spherical symme- 
try is dubious . The best study of Coma was performed by 
ICavazzi et al.l (l2009f ) , but unfortunately found only a rather 
speculative NFW profile of Af2oo = S.ll" x lO^*A/0 and 
'"200 = 1.81q gAfpc with no prior on the concentration pa- 
rameter (found to be C200 = S.Olj'g)- With a prior on the 
concentration parameter (set to be C200 = 3.5^|!| g) they 
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found M200 = 9.7l^:g X lO" a nd raoo = 2.2tni Mpc. Fi- 
nally, we added an estimate from lKubo et al.l |2007|) : M200 = 
27.0 ± 8.0 X 10", r2oo = 2.9 ± 0.3 Mpc and C200 = 3.8l"g^. 
The mass of the central galaxy would need to be 1.5, 1.0 and 

2.0 X IO^'^Mq respectively for the three cases. The K-band 
luminosity of the central galaxy in Coma is IO^'^Lq, so there 
is plenty galactic mass to supply what is needed. 

4.1.6 RX,J I4I6 

We included one fossil group RXJ 1416 in the sample since it 
has a very minor galaxy component (aside from the BCG). 
iKhosroshahi et al.l (|2006l ') fitted an NFW profile to X-ray 
data from Chandra and XMMN with parameters M200 = 

3.1 X 1O^*M0, r2oo = 1.2 Mpc and C200 = 11.2 ±4.5 and we 
used it to infer the necessary M/L of the BCG which has 
luminosity Lr = 7 x lO^^L©. The M/Lr need only be 1.2, 
whereas 5 is a typical value for an old stellar population in 
that band. So here again, as with several of the other massive 
systems, the full luminosity of the BCG is not required al- 
though the SN halo reaches the TG limit at the centre. This 
might be because the BCG is extended and only a fraction 
of the light is enclosed within 20 or 30 kpc (where the SN 
halo reaches the TG limit). Alternatively, if we used the full 
mass of the BCG, the SN density would simply reach the 
TG limit at a smaller radius. 



has a relatively young stellar population, although there is 
a age gra dient from ~ 12 Gyr in the centre to ~ 3 Gyr by 
~ 1 kpc (|Humphrev et al" II2OO6'). This leads to significant 
systematic uncertainty in constraining the M/Lk, which is 
given as 0.54 ± 0.11 and 0.35 ± 0.07 for the Salpeter and 
Kroupa IMF respectively. A recent re-evaluation with data 
of superior resolution (D. Buote private communication) 
puts the Kroupa value at 0.49 ± 0.18 (meaning the Salpeter 
value will be somewhat larger), but this still falls well short 
of the necessary galactic mass. This could be a serious 
problem if the low M/Lk could be confirmed, but for now 
it is a prediction of this model that the true M/Lk (when 
the correct IMF, age distribution and merger history are 
taken into account) will be close to 1.2. 

NGC 720 is not the only system that is very sen- 
sitive to the ob s ervati ons. NGC 1550 was studied by 
ICastaldello et al.1 l|2007l ') with both Chandra and XMM- 
Newton and both sets of temperature data are plot- 

ted in their Fig 3; Later it was observed again by 

iKawaharada et al] (|2009l ) using XMM-Newton, in a study 
which suggested there was evidence for a recent merger. If 
we use either set of XMM-Newton data points, the M/Lk 
would need to be at least 1.6, whereas if we only subscribe 
to the Chandra data, the M/Lk need only be 1.2. Follow up 
observations could provide a very strong test of the model, 
but Chandra's greater spatial resolution makes it the more 
reliable data set. 



4.L7 NGC 4125 and NGC 6482 

These two X-ray bright early type galaxies were taken from 
the sample ofiHumphrcv et al. (2006). Interestingly from the 
point of view of SNs in galaxies (which could disrupt the 
good fits to rotation curves i.e. the basis of MOND) their 
dynamical masses are comfortably covered by the BCG mass 
everywhere. From their density figures you can see the SNs 
are at the TG limit everywhere, but the SNs make virtually 
no impact on their mass profiles. In fact, if the SN density 
for NGC 4125 was 8 x lO'^Mopc"^ all the way to the centre 
from 10 kpc, the enclosed mass of SNs would be 3.4 x 10^ M© 
in comparison to roughly 1.0 x lO^^M© for the BCG (at that 
radius), which is a factor of 30. Therefore, there is no reason 
to believe that SNs will influence the internal dynamics of 
individual galaxies in a meaningful way. 

4.1.8 NGC 720 and NGG 1550 

These are the two most p roblematic groups. 
iBuote fc Canizared l|l994 Il996l ): iBuote et all (|2002l ') 
observed that the twisting of X-ray isophotes around the 
elhptical galaxy NGC 720 compared to the intrinsic ellip- 
ticity of the galaxy (which outweighs the gas by more than 
two orders of magnitude) could only be compatible with the 
presence of DM at least four times the galaxy mass (by the 
edge of the galaxy). This SN halo would have an ellipticity 
distinct to that of the galaxy and the gas would trace 
the superimposed potential, hence generating the twisting 
isophotes. The luminosity of NGC 720 is Lk = 1.7 x lO^L© 
and requires M/Lk = 1-2 which is the high end of the 
scale for the Kroupa IMF of an old stellar population. 
The SNs have reached the TG limit by 30 kpc, so cannot 
contribute any more mass. The problem is that NGC 720 



4.2 Sterile neutrino mass and the /^-function 

lAngusI (|2009l ) demonstrated that 11 eV/(? SNs were required 
to fit the first three acoustic peaks of the CMB and that 2 
species of 5.5 eV/c^ SNs were totally inadequate (as were 
three 2.2 eV/c? ordinary neutrinos). There was also very 
little freedom beyond 5 per cent of 11 eV/c^ , without intro- 
ducing unjustified free parameters. 

For further instruction to the importance of the SN 
mass to be 11 eV/c? , we calculated the M/Lk necessary to 
produce an equilibrium SN distribution in the group ESO- 
306 (which is by no means the most constraining) in two 
different scenarios. Firstly, trying a 9.5 eV/(? SN and sec- 
ondly keeping an 11 eV/cP' SN mass but using the standard 
^-function of MOND: where instead of using Eq|4]we try 
^.{g/ao) = , Whereas using the simple /i-function 

and an 11 eV/c^ SN requires M/Lk = 1, the 9.5 eV/c^ SN 
requires M/Lk ~ 1.5 and the standard /i-function would 
need 1.4. Given that these sorts of high M/Lk would be 
the rule, rather than the exception, they seem incompat- 
ible with stella r popu lation synthesis models like those of 
iBell fc de Jo"n3 (|200ll ). 

This evidence in fa vour of the sim ple ^-function is not 
an iso lated case. In Fam aev fc BinnevI (2005) and McGaughl 
(|2008l ). the simple function was prefered from a fit to 
the Milky Way's rotation curve (a high surface brightness 
galaxy) as was the case in the larg e sample of high surface 
brightness galaxies carried out by ISanders fc No ordcrmeeij 
(2007). This preference is also being found by an ongoing 
study of ultr a high resolution ro tation curves observed by 
the THINGS ([Walter et al.ll2005h collaboration (G. Gentile, 
private communication). The main problem with the sim- 
ple function is that it cannot be used all the way to the 
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strong gravity regime. In particular, it would produce too 
high a modification in the inner Solar System, which is ex- 
cluded, for example, by measures of the perihelion precession 
of Mercury. The solution is to use a /i-function that would 
rapidly interpolate between the simple and standard /j, for 
very large values of the gravitational acceleration, but this 
is not applicable here. 



4.3 Common central phase-space density 

It is an encouraging result that the 11 eV/c'^ SNs have 
unique and continuous equilibrium models for such a large 
dynamic range of cluster properties, but what is remarkable 
is that at the centre of every cluster, the TG limit is reached. 
For each given cluster, there is a specific, maximum density 
profile that can exist in equilibrium: a good example being 
NGC 5129 (in Fig [3}, which must have the maximum al- 
lowed density from 40 kpc for the SN halo to both be in 
hydrostatic equilibrium and provide the correct dynamical 
mass as measured by the ICM properties. 

Accordingly, it is the TG limit and the degenerate prop- 
erties of the SNs (since they are fermions) that sets the dy- 
namical properties of all clusters. No relation, even remotely 
like this exists if the cluster DM is cold or non-fermionic. 



4.4 Sterile neutrinos inside galaxies 

4-4- 1 Dynamics 

It is apparent from the examples of NGC 4125 and 
NGC 6482 (see H4.1.7|l that SNs wiU mo stly not aflfect the 
dynamics of galaxies. Gentile et al.l (|2008l ) have shown that a 
constant density of 10~^ Mqpc~^ was allowed in the MOND 
fits to the rotation curves of Ursa Major galaxies. This mag- 
nitude of SN density is typically found within 1 AIpc in the 
very massive clusters and within 100 — 300 kpc from the cen- 
tre of groups. At the centre of clusters, the situation is ob- 
viously different, but stable spiral galaxies are never present 
there since tides would rip them apart before they fall to the 
centre. Therefore, this will not disturb the MOND TuUy- 
Fisher relation since field spirals should be far from any SN 
halos. 

However, SNs are required to exist in the centres of some 
very massive ellipticals. This is clear from the right hand 
panels of the enclosed mass profiles of the clusters A 478 
(FigH) A 907 and others. In these figures one can see the 
mass of the BCG (the solid line) is smaller than the MOND 
enclosed mass. Thus, a considerable mass of SNs is required 
inside the limits of the stellar orbits, but always significantly 
less than the mass in stars. 



4.4-2 Lensing studies of individual galaxies in MOND 

In its original form, MOND should be able to explain all 
galaxy dynamics without dark matter. Although rotation 
curves have always yielded exc ellent results, the data for 
lensing studies (beginning with IZhao et al.l |2006| ) of indi- 
vidual galaxies h as not been as promising. For instance 
iTian et ahl (|2009| ) show that the weak-lensing of single, iso- 
lated galaxies is perfectly compatible with MOND up to 
a particular galaxy luminosity (from Lr — 0.1 — 8.0 x 
10^" L0). Thereafter, the lensing data implies the need for 



DM, which is exactly as we might expect if these more 
luminous galaxies are embedded in a low density but ex- 
tended SN halo (akin to those in the lHumphrev et al.ll2006l 
sample, see j 4^1.711. The same is true for the study by 



iFerreras et al. 



I 2008h whereby the integrated mass along 
the line of sight is entangled with the stellar mass. This 
makes lensing studies of individual early-type galaxies 
in MOND significantl y inferior to those using ICM, or 
even globular clusters (iRichtler et al.ll2008i). planetary neb- 



ulae iRomanowskv et al.ri2003l: iMilgrom fc SandersI |2003| : 



iDouglas et al. _ 20071 : 'Napolitano et al.1 200E ) and satellite 



galaxies ijKlvpin fc Prada 2009 : Angus et al 



l2008b[) . 



5 DISCUSSION AND CONCLUSION 

In this paper we have shown explicitly how to calculate the 
equilibrium configurations of neutrinos in MONDian galaxy 
clusters. The density of sterile neutrinos is fixed by the prop- 
erties of the ICM (or the observed lensing map), but deriva- 
tion of the sterile neutrino velocity dispersion allows a spe- 
cific density profile to exist in hydrostatic equilibrium. We 
have presented the detailed properties of 30 typical galaxy 
groups and clusters over a wide range of masses and tem- 
peratures (and even redshift vis-a-vis the bullet cluster, the 
lynx cluster and A 1689) and have shown that not only can 
we elucidate velocity dispersion profiles that allow the sterile 
neutrinos to exist in hydrostatic equilibrium, but also that 
the Tremaine-Gunn limit sets the central density. 

It would appear to be a very strong coincidence that 
by doing little more than fixing the mass of a sterile neu- 
trino to be 11 eV/c^, we can serendipitously, explain the 
formation of the acoustic peaks in the CMB and specify 
the exact properties of systems that require DM in MOND. 
In particular, regardless of cluster mass, the velocity dis- 
persion of sterile neutrinos necessary to impose hydrostatic 
equilibrium allows the density to reach its maximum (which 
is a function of both sterile neutrino particle mass, velocity 
dispersion and cluster mass indirectly) at the centre of the 
cluster. The only stipulation is whether these equilibrium 
configurations are stable, which should be the next check. 

Out of the 30 systems, there are two which need to be 
monitored. NGC 720 requires a K-band mass-to-light ratio 
of 1.2, which is possible for an old stellar population, but the 
current best population synthesis model suggests 0.49±0.18. 
Estimate s of stellar masses a re notoriously fraught with dif- 
ficulties l|Conrov et al.l [20091 ) . but this is a potential prob- 
lem. NGC 1550 (discussed alongside NGC 720 in H4.1.8|l can 
only reproduce the temperature observed from Chandra; if 
the XMM-Newton data are used, the K-band mass-to-light 
ratio of the central galaxy is too high. From this study, it 
is clear that the strongest tests of the 11 eV/c^ sterile neu- 
trino model comes not from the rich clusters, but rather 
from smaller groups or individual galaxies with bright X- 
ray halos. 

In a similar sense to how the NFW density profile de- 
duced from N-body simulations of CDM structure forma- 
tion have been shown to be inadequ ate descriptions of some 
galaxy and galaxy cluster DM halo s llBroadhurst fc Barkanal 
l2008l : Ide Blok fc McGaughl Il998l : ICentile et all |2004 ) we 
suggest that our highly regular density and velocity disper- 
sion profiles shown in Figs IH6I must be used to judge the 
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consistency of simulated structure growth in MOND cosmo- 
logical simulations. Presumably, it is guaranteed that the 
objects that condense out of the background will reach the 
Tremaine-Gunn limit at the centre (which makes it so al- 
luring that all the clusters studied here do), but it is not 
certain that the phase space densities will fall with radius 
in the manner shown here, nor that the chemical potential 
will have the correct values. 

Given that we fit the CMB and have these interest- 
ing results for clusters of galaxies, another crucial test of 
this model (after checking stability) will be to see if MOND 
N-body cosmological simulations, with 11 eV/c^ sterile neu- 
trinos can form structures resembling those shown here and 
match th e linear matter power spectrum. It is worth pointing 
out that ISkordis et al.l (|2006l ^ computed the matter power 
spectrum with MOND-like gravity and three 2.75 eV/c^ 
active neutrinos (the ao used was 3.5 x larger than the 
one used typically to fit rotations curves and used here) 
and showed the extra matter density from the neutrinos 
(f2„ — 0.17 compared to Sit ~ 0.05) coupled with the 
MOND gravity could come relatively close to providing a 
good match to the matter p ower spectrum measured by the 
SDSS (jTegmark et al.ll2004 ). The higher matter density pro- 
vided by the SNs (O^^ = 0.225 and Qt = 0.047) traded off 
with the smaller, more standard ao, could very well provide 
a superior fit to the matter power spectrum. To resolve this 
we need cosmologic al numerical simulati ons, which are still 
in their infancy (see iLlinares et ahlbOOSl ). 

Relativistic MOND theories begi nning with T eVeS 

(jBek cnstcin 2004 , also see the review bv lSkordidlioogl ) and 
BSTV ( Sanders! 12^0051 ). which le d to new ideas Uke Gener - 
alised Einstein- Aether (see e.g., IZlosnik et al] |2007| . 120081 ) 
are still highly complex and only address the galactic dark 
matter problem (not the cluster, or cosmological dark mat- 
ter problem; nor the dark energy problem). Therefore, one 
caveat we would add to our conclusions is that this cosmo- 
logical model still requires dark energy in the same coinci- 
dental amount as A. CDM. H o wever, as elaborated upon in 
iMilgromI (jl999t ) and lMilgroiiil (j2008ah . MOND and dark en- 
ergy must be two sides of the same coin that leads seamlessly 
to 2- Rao ~ c(A/3)^^^. There is some progress in this direc- 
tion dFiizfa fc Aiimill2007l; iBruneton et"aLll2009l: ITi fc Zhaol 
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Blanchet fc Le Tiediioogl ) but it has not yet been con- 



vincingly or efficiently shown. 



Table 1. Here we list the K-band luminosities of our BCGs along 
with the minimum K-band mass-to-light ratio required to fit the 
dynamical mass. The clusters with a (*) lack information about 
the BCG, so in the luminosity column, we have entered the re- 
quired BCG luminosity with unity M/ L. The clusters with a lumi- 
nosity followed by (V),(B) or (R) have their luminosity measured 
in that band, and not the K-band. The Coma cluster has 3 sepa- 
rate mass profiles, so the 3 M/Ls refer to the profiles in the order 
they are taken in Fig|5] The clusters are separated into samples; 
the top set are from the Gastaldello et al. (2007) sample; then the 
Vikhlinin et al. (2006) sample; Humphrey et al. 2006; Zappacosta 
et al. (2006); miscellaneous NEW fits; the bullet cluster. 
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Figure 1. As per Figs|2]to|6] each row refers to a specific cluster, whereas the columns are properties of that cluster. In the left hand 
panels we plot density against radius for sterile neutrinos (solid), intracluster medium (dotted) and Tremaine-Gunn limit (dashed). In 
the central panels we plot velocity dispersion against radius for sterile neutrinos (solid) and intracluster medium (dotted). In the right 
hand panels we plot enclosed mass against radius for the Newtonian Dynamical mass (dashed), MOND dynamical mass (dot-dashed), 
intracluster medium (dotted) and the total mass of the brightest cluster galaxy (solid). Each cluster's designation is inscribed on the left 
hand panel. 
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Figure 2. As per Fig[T] 
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Figure 3. As per Fig[T] 
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Figure 4. As per Fig [T] except certain lines corresponding to the intra-cluster medium are absent because the sterile neutrino densities 
were deduced from weak or strong lensing analyses. 
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Figure 5. As per Figs I H4l except that tliey are for A 1689 (Halkola et al. 2006), tiie Coma cluster, tlie two clusters tiiat comprise the 
Lynx cluster (Jee et al. 2006) and the fossil group RXJ 1416 (Khosroshahi et al. 2006), from which the masses are taken in NFW form. For 
the Coma cluster we use 3 different measurements for the NFW mass profile as discussed in detail in i|4.1.1l To avoid crowding the Coma 
density plot, we only plot the sterile neutrino density for NFW profile 1 (solid), 2 (dotted) and 3 (dashed) and not the Tremaine-Gunn 
limit, which is instead marked by an asterisk, plus sign and diamond respectively where the sterile neutrino density reaches it. 
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Figure 6. As per Figs lHSI cxcept that they are for the two clusters that comprise the bullet cluster llClowe et al.ll2006l) . from which the 
masses of the two clusters are taken in NFW form. This case is discussed in detail in ^4.1.11 
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